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C^ ' Abstract 



a^ 



Variation of coupling constants of integrable system can be considered as canonical transfor- 
mation or, infinitesimally, a Hamiltonian flow in the space of such systems. Any function T(p, q) 

CO . generates a one-parametric family of integrable systems in vicinity of a single system: this gives an 

^ ' idea of how many integrable systems there are in the space of coupling constants. Inverse flow is 

'^ , generated by a dual "Hamiltonian" , T(p, q) associated with the dual integrable system. In vicinity 

of a self-dual point the duality transformation just interchanges momenta and coordinates in such 

r — , a "Hamiltonian": T{j),q) — T{q,p). For integrable system with several coupling constants the 

f^ ' corresponding "Hamiltonians" Ti {p, q) satisfy Whitham equations and after quantization (of the 

original system) become operators satisfying the zero-curvature condition in the space of coupling 
constants: 

- Ta{p, q), Tb(p, q) 



o 
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Some explicit formulas are given for harmonic oscillator and for Calogero-Ruijsenaars-Dell system. 



-J. 1 Introduction 

X 

C^ I The study of non-perturbative quantum phenomena revealed the real role of integrable systems in 
physics: r-functions of integrable hierarchies appear to describe partition functions of quantum the- 
ories, with time- variables identified with coupling constants [ffl, 0]. Broadening understanding of this 
relation, together with the growing interest to peculiarities of essential (non-perturbative) quantum 
physics, stimulates a new attention to the basics of integr ability theory, which brings one back to the 
very foundations of classical Hamiltonian mechanics that were not sufficiently investigated in the past. 
The present paper is devoted to one of such basic subjects, to the description of integrable systems in 
terms of the Hamiltonian evolution in the space of coupling constants. This approach can help one to 
describe the entire variety of integrable systems and to reveal the meaning of peculiar relations, like 
dualities, between different such systems. 

2 Hamiltonian flow in the space of integrable systems 
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2.1 Variation of couplings as canonical transform — peculiar feature of integrable 
systems 

Integrable system with A^ coordinates qi and N momenta pt is characterized by existence of A^ Poisson- 

commuting Hamiltonians Hi{p,q), {Hi,Hj} = -^-gf- M'~U^ ~ ^- ^^^ such a system one can 

consider a canonical transformation, treating these Hamiltonians as new momenta-like or coordinate- 
like variables. This transformation (of which the infinitesimal version is a certain Hamiltonian flow) 
will be the subject of the present paper. 

To make the problem precise, let us consider a one-parametric family of integrable models, parame- 
terized by a single coupling constant g such that the model is free when g = 0. This means that at g = 
the Hamiltonians H^ (p) = Hi{p,q\g = 0) are functions only of momenta p, though, for conventional 
choices of Hamiltonians in particular applications, these functions can be non-trivial. The typical 

examples are: Hj^ (p) = J2i=oPi ^^'^ ^k \P) ~ J2r.\i\=k'n.i£i ^^' fo^' p-rational and p-trigonometric 
models respectively and some elliptic functions of p for their elliptic generalizations. 

The adequate definition of the new canonical variables Pg = P{p, q\g) and Qg = Q{p, q\g)l\ is @ 

Hl'\p)=H,{p,g]g) 

Hf\q)=H,{Q,P\g) 

where Hi{p,q\g) define the Hamiltonians of the dual integrable system (for other discussions of this 
duality see [Q] and |p). Note that in (||) they depend on the "dressed" variables P and Q and, moreover, 
P and Q are interchanged. The shape of the dual Hamiltonians is dictated by the requirement that 
the new variables P and Q are canonical, i.e. 



y^ dPi A dQi = ^ dpi A dqi 



and the Poisson brackets are 



/ d d d d \ ^^ f d d d d \ 

2.2 Infinitesimal version of canonical transforms — the Hamiltonian flows 

Relations (|l]) define Pg as functions of p, q and the coupling constant g. Of special interest and 
importance is the infinitesimal version of this canonical transformation, considered as a Hamiltonian 
flow along the g direction in the space of coupling constants. Such transformation is generated by a 
new Hamiltonian T(j), q; g), according to the rule: 

^ = {T(4<?,|rf,0,} = || 



^ In what follows we often omit the label g, implying that the capital letters P and Q always denote the dressed 
momenta and coordinates Pg and Qg. 



This new Hamiltonian does not commute with the old ones, but it converts them into a new set of 
commuting Hamiltonians. 

Moreover, instead of considering a pre-given family of integrable systems, one can use any function 
T{p,q\g) to define a whole one-parametric family of integrable systems, though explicit construction 
of the corresponding Poisson-commuting Hamiltonians is rarely possible. Further, a multi-parametric 
family of integrable systems is similarly generated by any collection of functions Ta{p,q\g) satisfying 
the compatibility (Whitham) equations:^ 

dTh dTa , , 

dga dgb 

The coupling constant variation for the dual integrable system is governed by the dual Hamiltonian 
f{p,q\gy. for 

Hf\P)=H:{p,q\g) (6) 

(and Hf^\q) = Hi{Q,P\g)) we have 

(7) 



Relation between T{p,q\g) and T{p,q\g) becomes especially simple at any self-dual point (where 
Hi{p, q Igso) = Hi{p, q \gsD))- 

T{p, q \gsD) = T{q,p \gsD) (8) 

As immediate consequence, any symmetric function T{p,q\g) = T{q,p\g) defines a one-parametric 
family of self-dual integrable systems. 

Note that it is not that immediate calculation to restore the Hamiltonian starting from arbitrarily 
given T{p,q\g). Say, one can build a perturbative series in the coupling constant g. At the self-dual 
point for the system with one degree of freedom one takes the expansion T{p, q\g) = J2i=i Ti{p, q)g^~^, 
P = H{p,q) =p+J2i=i^i{P,Q)9' etc and obtains that ^iip,q) = dgTi{p,q), ^2ip,q) = dqT2{p,q) + .... 
The result ultimately is that the symmetric part of dp^i{p,q) = dgpTi{p,q)+ is given independently 
at any order, while the antisymmetric parts are fixed by $j's at lower orders: 

^ As mentioned in the abstract to this paper, after quantization of the original system these Ta{p, q\{gb}) become a 
family of operators Ta{p, <l\{gb}), satisfying the zero-curvature condition in the space of coupling constants: 



dg 


Yyi'g,l^g\g),I'gJ - 


dQg 


dQg 
dg 


= |f(P^,cJ|5),?,} = 


df 

dPg 



dga dgb 



= 



This can serve as an explanation of emergency of classical integrability (the zero-curvature equations) in the study of 
quantum integrable systems. 



<^[{p,q)-^[{q,p)=0, 

^'3ip,q)-^'siq,p) = ^i{p,q)^2{p,q)+mp,q)^l{p,q) - h'l'{p,q)^'iip,q)- (9) 

-^i{P,q)<^'i{p,q)<^i{p,q) - {p ^ q), 



where all the derivatives are taken w.r.t. to the first variable. A trivial solution to these equations 
(|9|) is ^i{p,q) = p"*+^gi"» with arbitrary {tij}. This implies that the Hamiltonian H{p,q) = pf{pq) 
is self-dual with arbitrary function /. Indeed, one easily checks this is the case (note that for such a 
Hamiltonian pq = PQ). We discuss later in detail the rational Calogero Hamiltonian that gets to this 
class of Hamiltonians. 

2.3 Generating functions of canonical transforms 

Along with the Hamiltonians T{P,Q\g) one can also consider the generating functions of canonical 
transformations in question, like S{Q, q\g) or its Legendre transform F{P, q\g) = PQ — S{Q, q\g), such 
that 

-^ ds ^ ds 

-P=^, P=T7^ (10) 

dQ dq ^ ^ 



and 



Of course, for canonical transformation 



-f 



^ + ^ = (12) 

as implied by (|^), but one should be more careful when drawing similar conclusion from (11): the sec- 
ond g derivatives satisfy eq.(g), because g'-derivative is taken at constant Pg and Qg, which themselves 
depend on g. 

Similarly, 

-^ dF ^ dF 

At self-dual points F{P,q\gsD) = F{q,P\gsD) is a symmetric function. 
Note that 

dF 

(with p and Q constant), i.e. T is, in a sense, a more invariant quantity than S and F, which does 
not depend on the choice of independent variables. 



Another way to see it is to consider a flow from an integrable system with couphng constant gi 
to the same system with coupling constant 52- Then, T depends only on §2, but not on gi, while the 
generating functions depend on both g^ and 52- 

2.4 Generating functions, evolution operator and eigenfunctions 

Basically, if reexpressed in terms of Q and q (instead of P and Q) the Hamiltonian T = dS/dg. 
Therefore, e can be considered as a kind of an evolution operator (kernel) in the space of coupling 
constants, which performs a canonical transformation from the free system to the integrable one. After 
quantization it can be symbolically represented as 

e'HQ,Q\g) = ©AlV^f (Q))cx(V'A(g1l (15) 

where I'^a) and 1-0^ ) ^^^ eigenfunctions of the system and cr are some coefficients, depending on the 
spectral parameter A. The evolution operator satisfies the quantum version of eq.(||), 

i/(°) {d/dQ) = e-'^^^^^3)H{d/dq, glg)e^^(^'«l3) (16) 

Since eigenfunctions \'il)\ (Q)) of a free system, are just exponents of the spectral parameter A, the 
dressed eigenfunctions {ipxiq)) are basically Fourier transforms of the evolution operator e^^'^''^': 

Mq)^ fe'^^^'^'^e'^'^dQ (17) 

Similarly, 

Mq) ~ e*^(^'«)5(A - P)dP ~ e*^(^'^) (18) 

Note that there is a freedom in solutions of eq.(2) to shift {Qi} by an arbitrary function of {Pi}- This 
shift is quite complicated in terms of the generating function S, but in F it is just an addition of the 
term depending only on {Pi}. In the quantum case, this ambiguity in the definition of F is just a 
matter of normalization of the eigenfunction. 

One can also consider the quantum counterpart of T{p, q\g) which is a Hamiltonian that gives rise 
to a Schrodinger equation w.r.t. the coupling constant 

'^=^n (19) 

Therefore, the wavefunction ■0 can be also realized as a path integral over the phase space variables 
P{9),Q{9)- 

The generating functions S{Q,q), F{P,q) and T{p,q) satisfy the quasiclassical versions of these 
relations (when multiple derivatives of S, F and T are neglected). Exact definition of quantum 
evolution operators, including effects of discrete spectra and identification of the spectra at different 
values of coupling constants, i.e. precise definition of the spectral parameter A in such a way that 
eq.(ll5|) is diagonal in A, will be considered elsewhere. 



3 Some conceptual questions for further examination 

The previous section explained possible needs to study integrable systems from the non-conventional 
perspective: in terms of the auxiliary Hamiltonians T{p,q\g), describing the associated Hamiltonian 
flows in the space of coupling constants. Some immediate problems were also mentioned, like 

(i) criteria, distinguishing the Hamiltonians T associated with conventional integrable systems at 
the entire space of such Hamiltonians; 

(ii) structure of algebra, generated by all the Hamiltonians Hi and T^, in which the ffj's generate 
a commutative (Cartan) subalgebra; 

(iii) duality relation between T and T; 

(iv) relation between integrable systems and flat connections over the space of coupling constants 
encoded in (|5|). 

One can also add to this list: 

(v) identification of the flows and Hamiltonians Ta in terms of Seiberg-Witten theory, |^, |2|; 

(vi) treatment of T's in terms of the theta- function formalism of refs.|^, ^; 

(vii) the "gauge-equivalence" of the Hamiltonians Ta, corresponding to different choices of the 
constant-p, q sections over the moduli space (i.e. with different quantities kept constant when the 
couplings ga are varied); 

(viii) the role of the distinguished "flat-modulus" a Q, for which the WDVV equations Q hold; 

(ix) further development of the notion of generalized r-function, making use of the evolution 
operator e"^, introduced in the previous section (see also the concluding section]^ below). 

Since the subject is relatively new, we do not think it is time now to get into further discussions of 
these (and many other) problems. Instead in the remaining part of this paper we consider some simple 
examples: the Hamiltonians T{p,q;g), associated with the flows in the coupling constant space for 
some well-known integrable systems like harmonic oscillator and Calogero-Ruijsenaars-Double-elliptic 
(DELL)0 family. 

4 Particular examples of Hamiltonian flows 

4.1 Harmonic oscillator 

Perhaps surprisingly, even for this simple model one obtains rather sophisticated formulas. 

For harmonic oscillator H{p, q\uj) = ^{p"^ + uj'^q'^). Let the frequency lo play the role of the coupling 
constant so that H^^'{p) = ^p'^. Then 



P = Jp^ + uP'q^ , 
Q = arctan — 



uj p 

Note that in this case the w-evolution possesses "the conservation law" 



^ In fact, the treatment of the DELL case is very sensitive to the choices mentioned in (vii) above, and we postpone 
serious discussion of the DELL case to another occasion, in order not to draw attention away from the simple ideas of 
the present paper. 
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The inverse transformation looks like 



p = P cos 
P 



Loq 
P 



ujQ 
ujQ 



(21) 



(22) 



— sm ■ „ 
w P 



i.e. the Hamiltonian of the dual flow is 



H{p,q\to) 



q ujp 

— sm —zr 
uj a 



(23) 



(Note that with this definition H^^' (p) = p and also note that p, g-duality does not respect conventional 
dimensions of p and q so that it should not be a surprise that ujq/p in H got substituted by cop/q in 
H. Of course, p and q are nothing but Q^^ and P^j.) 

The generator of a;-evolution such that 



dP 
dP 



ujq 



P 



p,q=const 



„2 



■ sm 



uj(4 



dT 



p- 



p,q=const 
P 



<JJ^\/P^ _|_ ^2g,2 



ujq 

p 



arctan 



bjq 
p 



IS 



ijj 



T 



2 UjQ ( OjQ LoQ 

2 COS — (tan-^ — 



dT 

ap' 



P2 2loQ PQ pq-PQ 

— t: sm — — = 

4u;2 P 2uj 2uj 



For the dual system we similarly have 




ijjp ujp 

— ^ cos — ^ 

q q 



sm ■ 



arctan ■ 



Q . 



ujp 
Q 

f 

Q 



and 



(24) 



(25) 



(26) 



-LOPQ + (Q2 + ^2p2) j^^p^j^^ a£ 
2ZJ2 

pg - PQ PQ - pq 



2oj 



2uj 



-T 



(27) 



In order to obtain the generating function S{<Q, g|c<j) one needs first to reexpress p and P through q 
and Q. Unfortunately, the formulas are transcendental. Instead, it is easy to express p and Q through 
q and P, this involves only elementary functions 



The corresponding generating function 



p = \J P"^ — uj'^q^, 

P u:q (28) 

Q = — arctan ■ 



F(P, q\uj) = / pdg = / Y P^ — oj'^q^dq 



so that the relation (|l^ ) in quasiclassical approximation is 

AV2w 



,2„2 



(29) 



\ A V A2 y (30) 

~ e-^'?'/2//e^(^g) 
where 2w ~ '^ ~'~ 2 ^^*^ ^'^^ Hermite polynomials, satisfying 

{-dl + x2)e-"'/2^e,(x) = {2u + l)e-"'/2^e,(x) (31) 

are given by inverse Laplace transform of the generating function e~^ l^+th 

He.{x)^j^,e-''/'^'^ (32) 

Taking this integral in the saddle point approximation, one gets 

He,[x) r.(x- ^x^-2{u + l)\ "^^ e^V2+x/2Vx^-2(.+i) (33) 

which, at large v, is in accordance with eg. (|30|) . 



4.2 Generic Schrodinger-like systems 

For a single particle any dynamics is integrable, and any Hamiltonian is canonically equivalent to the 
free one. Therefore, it may make sense to look at the corollaries of our formalism in such context. This 
can help to see the difference between models which have and have not simple integrable generalizations 
to the multiparticle case. 

For H = yjp^ + g^V{q) the Hamiltonian T can be found in the form of a perturbation expansion: 

npM9)=9^^ + g'^ + 0{g') (34) 

p p-^ 



where 



^...„, ?m--\T.^ 



There is no way to integrate explicitly these equations in the general case. However, for harmonic 
oscillator exact answer exists for the entire expansion, eq.(^5|). Similarly, exact answers are available 
for potentials of the rational and trigonometric Calogero models (see below). However, for elliptic 
Calogero model with V{q) ~ sn~'^{q\k) eq.(35) implies that the generator T{p, q\g) involves the integral 
elhptic tangent (i.e. the Weierstrass zeta- function) , Tn{q\k) = J gj^-iLf.) ~ C('?l''")-Q Higher orders of 
the g expansion contain multiple integrals of Tn(^). This reflects a more sophisticated group structure 
of elliptic models. 

The simplest way to derive formulas (|3^-(|35[) is to construct the Hamiltonian q = H{P, Q) dual 
to i/p-^ + g''^V{q). The canonicity condition then gives 



dH{P,Q) _ p{P,Q) 
dQ P 

Therefore, the dual Hamiltonian can be obtained via solving the equation 



(36) 



'^-I'TT^wmm '^'' 



w.r.t. q. Now taking the derivative of P = yjp'^ + g^V{q) w.r.t. g, one obtains 

-g = |v(,(P.Q)) (38) 

where q{P,Q) is given by (|37|). One easily obtains the expansion of these two formulas in powers of 
g, in particular, reproducing (p3)-(p5h. 



4.3 SU{2) rational and trigonometric Calogero-Ruijsenaars models 
4.3.1 The rational-rational case (rational Calogero model) 

In this case the single Hamiltonian is H = ^ (p"^ -\- ^j, thus H^{p) = ^p'^, and 



q^ 



(39) 



The rational Calogero model is self-dual. 



* However, the conservation laws of the (;-flow involve the ordinary (algebraic) elliptic tangents, tn[(^\k) = 
sn(^|fc)/cn(^|fc): say, in the 5*17(2) elliptic-rational model (the dual of the elliptic Calogero system) pi 

cn{P\k) = a{q)cn{f3{q)p\j{q)), 



sn2(P|fc)' 



«'(?) = 1-4' f3\q) = k" + k'a'{q), l{q) ^ ^, k" = 1 - k\ 

the counterpart of (Ell) is 

Q tn(P\k) = q tn(/3p|7). 

The difference between the algebraic and integral tangents is that ■^tn{^\k) = " .J, . while ■^Tn{^\k) = — tttttt 



One can express the dressed variables P = Pg,Q = Qg through the bare ones, p, q 

P = -Jp2q2 _ g2^ 

q V 

2 (40) 

Q 



vrV-ff^ 



The inverse transform is 



Note that, similarly to the case of harmonic oscillator, the g-flow possesses a conservation law 

PQ = pq (42) 

The flow is generated by the Hamiltonian 

rrf ^ li PQ-9 1, PQ-9 ..o^ 

T{p, q;g) = ^ log — -— = - log 43) 

2 pq + g 2 PQ + g 

This is obviously a symmetric function of p and q, thus, T£){p,q) = T[p,q) as it should be for the 
self-dual system. 

In order to obtain the evolution operator e , one should rewrite T in terms of the g, Q- variables. 
The p, P- variables are then given by 



,2_ 9' 



p 



Q^-q^' 

g'Q' 



(44) 



g2(Q2_^2) 

The symmetricity (integrability) condition (|l^) is true 



dq VQ' - q^ dQ q^Q^ - q^ 
and one obtains 



d I d Q 



S = g arccosh— = g\og\ = 7^ log , 

q \ q J ^ Q + VQ^ - q^ 

dSiQ,q\g) rj., i^ . I ^ 
Q = T(p{Q, q),q\g) 

In this case, S" is a simple linear function of g, and dS/dg = S/g. 
Similarly, 



(46) 
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^ ^ ^ 2 ^ypVTF + 5 (47) 
^ = T{p{P, q),q\g) 

One can also consider the transformation from the system given at {gi,Qi,Pi) to that given at 
((72, Q2, P2)- The corresponding generating function 



^,^ ^ , , 9l^^Ql-Ql + ^^9lQl-9iQl , 92^Ql-Ql + ^^glQl- glQl 

S{Ql,Q2,9l,92) =51 log ^^ -pr^ h 52 log -^ 

(48) 
is the function of both gi and g2, while 

55^1_^Q^-^ (49) 

is the function of only (72 • 
In accordance with (|l^) 

V'A(g)~ /e*^(«'«l9)e^Q^dQ~gJ,g(ag) (50) 



where Jg(2;) is the Bessel function. The exact quantum wavefunction, solving the equation (—9^ + 
^)'ipx{x) = \^ijj\{x), is ^/xJy{i\x), v^ = —g^ + 1/4. It coincides with ( |50D in quasiclassical approxi- 
mation. 

For other members of the Calogero-Ruijsenaars family we just list some relevant formulas, parallel 
to those in the rational-rational case. 

4.3.2 The rational-trigonometric case (trigonometric Calogero model) 

p2^/ 9^^ 

sinh q ' 

cosh^ q = cosh^ Q M ~ "52 ' 

P tanh Q = p tanh q, 

^ 1, PtanhQ — 5 1, p tanh q — g 

T = — log = — loET ; 

2 PtanhQ + 5 2 ptanhg + c/ (51) 

p2 _ 9'^ cosh^ Q 



sinh Q — sinh g 



2 g^ sinh Q cosh (7 

sinh (/(sinh Q — sinh q) 
cosh Q (9 sinh Q cosh g 



^9 . /«inT.2n_«lnT.2„ ^Q 



sinh Q ~ sinh g ^ sinh gy sinh Q — sinh g 
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In this case S is still a simple linear function of g 



S = g arcsinh — ) (52) 

V Sinn q 



4.3.3 The trigonometric-rational case (rational Ruijsenaars model) 

2 , 2 / ^i'^^ ^ \ 

cosh P = cosh pi K — , 

V 9 / 

sinh P 
QtanhP = qtanhp, 

1 QtanhP — tanhe 1, gtanhp — tanhe 
T = — log — — = — log ; 

2 QtanhP + tanhe 2 gtanhp + tanhe (53) 

sinh e 



sinh P 



Q^-q"' 



2 sinh eQ 

smh p 



{Q"^ — q^){q^ — smh e) ' 

^ ■ u 9 d gQ 

—- arcsinh — , „ + — — - arcsinh- 



dq ^Q2 _ q2 dQ 7(Q2-g2)(^2_sinh2e) 

In this case, S is no longer a simple linear function of the coupling constant e. 

4.3.4 The trigonometric-trigonometric case (trigonometric Ruijsenaars model) 

2 T 2 / sinh e\ 

cosh P = cosh pi 2~ ' 

\ sinh q I 

T 2 , 2 ^ / 1 sinh^ e \ 

cosh g = cosh Q \ 1 2 — ' 

Y sinh PJ 

tanhPtanhQ = tanhptanhg', 

1 tanh P tanh Q — tanh e 1 tanhptanhg — tanhe 
T ^ — log ^ — log ' 

2 tanh P tanh Q + tanh e 2 tanhptanhg + tanhe' (54) 

2 n sinh ecosh Q 
sinh P — 



2 „' 



sinh Q — sinh g 

2 sinh e sinh Q cosh q 

smh p — 



(sinh Q — sinh (7) (sinh g — sinh e) ' 



d sinh e cosh Q d sinh e sinh Q cosh q 
— arcsinh — - + — — arcsinh — 

^^ Jsinh^ Q - sinh^ g ^^ ^/(sinh^ Q - sinh^ g)(sinh2 g - sinh^ e) 
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4.4 SU{N) rational Calogero model 

To give an illustration of what happens in the niultiparticle case we start here with the simplest case 
of the rational Calogero model. More profound examples will be considered elsewhere. 

The Hamiltonians Hi{p, q) = TrU are defined as traces of powers of the Lax matrix Lij = piSij + 
g—. Accordingly Hf{p) = J2jP)- The Hamiltonian flow within the family of rational Calogero models 
(along the g direction) is governed by 

T{p, q;9) = -5E rV + ^(^') (55) 

jj Ptj'dij 

Technically, higher order terms can be done in the following way. As the first step, one needs to find 
dual Hamiltonians. Note that dual Hamiltonians technically each time emerge as necessary ingredient 
in constructing T. Moreover, with these intermediate objects, constructing T becomes well-defined 
although often complicated procedure and reduces to finding solutions of some partial differential equa- 
tions. Indeed, first of all the dual Hamiltonians emerge as solutions to partial differential canonicity 
equations and then T is obtained by integration the following set of equations: 

^ (dHk dT _ dHk dT\ dHk _^ 

These equations just follow from the fact that the dual Hamiltonians Hk are nothing but independent 
of the coupling constant g coordinates {qi\. 

Now, having H^ calculated, one can easily construct a recursive procedure of calculating T as 
series in (7: T = "^g ~^^Tk. Indeed, note that, in the leading order in g, H^ are just simple functions 
of Qi, say, just Qk which do not depend on g. This means that plugging the power expansion of T 
into (^6[) gives a system of A^ linear algebraic equations for the derivatives D^ = dTi/dPi. Solving 



these equations, one should just integrate D\ , finally obtaining Ti. Then one can use this calculated 
Ti in order to find T2 etc. Note that each time there emerge non-trivial integrability conditions 

dDi/dPj = dDj/dPi (57) 

Thus, the most difficult part of the calculation is to find dual Hamiltonians. In the rational Calogero 
system, when the system is self-dual, the calculation is very immediate. Similarly, the calculation can 
be easily done for any member of the Calogero-Ruijsenaars-Dell family. 

5 Conclusion 

A peculiar feature of integrable systems is that interaction is introduced by a canonical transformation 
of a free system. In other words, the corresponding flows in the space of coupling constants are Hamilto- 
nian ones. This observation should provide new insights into the structure of exact (non-perturbative) 
partition functions, which are always (generalized) r- functions of some integrable hierarchies ||l|, ||. 

Let us remind that the partition function, arising after functional integration over fields. 
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J Dipexp lA{ipo + ip) + i^t{k}V{k}{^o + ^) 



{k} 

(58) 




JS\ 



^{t{k}\fo) = (G((^o)-Pexp i}_^ / H{k}dt{k} ) = {Ge 



is essentially unsensitive to the choice of original action A{(p), and the equivalence classes (of ac- 
tions) can be represented by integrable systems, i.e. by the systems where the partition function is 
a matrix element of a group element G of some group |^]. Our evolution operator in the space of 
coupling constants, e = -Pexp ( z^r^j / f/'r^idtr^j j, which performs a canonical transformation of 
integrable system into the free one, as discussed in s.g should also be a group element. The moduli c^O) 
parameterizing the boundary conditions and vacuum states of the original theory are related to the 
(eigen) values of the Hamiltonians -f^{A:} (i-e., in our simple examples, the dressed momenta P). The 
couplings (perturbations) t x^i (i.e. our g) become time-variables of the r-function. The fact that the 
space of relevant Hamiltonians (the moduli space) is an orbit of the canonical transformation group, 
which can be alternatively considered as changing the couplings can further clarifies the relation be- 
tween the times- and moduli-dependencies of the partition functions 0] and bring us closer to the 
long-awaited group-theoretical treatment of entire r-function, its time-dependencies being associated 
with the commutative (Cartanian) part of the algebra (generated by the Hamiltonians Hi), while 
moduli-dependencies - with its non-Cartanian part (generated by Ta and their commutators among 
themselves and with Hi). See refs.[10| for related approaches to this problem. 
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